Different Hagedorn temperatures for mesons and baryons from experimental mass 
spectra, compound hadrons, and combinatorial saturation 
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We analyze the light-flavor particle mass spectra and show that in the region up to ~ l.SGeV the 
Hagedorn temperature for baryons is about 30% smaller than for mesons, reflecting the fact that 
the number of baryon states grows more rapidly with the mass. We also show that the spectra are 
well reproduced in a model where hadrons are compound objects of quanta, whose available number 
increases with mass. The rapid growth of number of hadronic states is a combinatorial effect. We 
also point out that an upper limit on the excitation energy of these quanta results in a maximum 
number of hadron states that can be formed. According to this combinatorial saturation, no more 
light-flavor hadron resonances exist above a certain mass. 

14.20.-c, 14.40.-n, 12.40Yx, 12.40Nn 



In 1965 Hagedorn [y postulated that for large masses 
m the spectrum of hadrons grows exponentially, p{m) ^ 
exp {m/TH), where Th, the Hagedorn temperature, is a 
scale parameter. The hypothesis was based on the ob- 
servation that at some point a further increase of energy 
in pp and pp collisions no longer raises the temperature 
of the formed fireball, but results in more and more par- 
ticles being produced. Thus, there is a maximum tem- 
perature that a hadronic system can achieve. The Sta- 
tistical Bootstrap Model predicted that asymptotic 
behavior, namely p(m) ~ cm°' exp (m/TH), where a is 
a negative power (a < —5/2 [||, or a < — 3 |^). The 
parameters of the asymptotic spectrum are to be deter- 
mined by comparing to the experimental mass spectra. 
The fits of the sixties had a rather poor spectrum to their 
disposal, sufficiently dense only in the range of masses up 
to IGeV, and sparse above. Still, Hagedorn and Ranft 
[01 were able to fit the function 
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with mo = 0.5GeV, and got Th = 160MeV. Frautschi |] 
obtained similar values for the Hagedorn temperature. 
Ever since it was believed that there is one universal scale 
in the asymptotic spectrum of hadrons, and its value is 
about 160MeV. Meanwhile, the particle tables became 
more complete and longer, with a total of 3182 light- 
flavor states (counting mesons, baryons, and antibaryons 
with their spin and isospin degeneracies) compared to 
1432 used in Ref . [Q . This allows for a much more strin- 
gent verification of the Hagedorn hypothesis. We have 
done this, with surprising results. 

We have used the 1998 edition of the Particle Data 
Group tables |^ . Rather then comparing the experimen- 
tal and theoretical spectra, we compare their cumulants: 
the number of states with the energy less than to. 
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where is the step function, is the mass of the reso- 
nance, and gi is its degeneracy. Obviously, dN (m) / dm = 
p (to) . Although the information content is the same 
whether one uses p{m) directly, or N(m), the use of cu- 
mulants is advantageous, since that way one avoids the 
nuisance of producing histograms of pcxp- Our results 
are presented in Fig. 1, where we show the logiQ N{m) 
separately for mesons and baryons. We can clearly sea 
that up to TO ~ 1.8GeV the data line-up along (almost) 
straight lines. The solid lines are the least-square fits 
to the data according to Eq. (|l]) over the range up to 
1.8GeV, skipping the lightest particle in the set. The 
behavior of Fig.l is compatible with the Hagedorn hy- 
pothesis. 
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FIG. 1. Cumulants of meson and baryon spectra, Hagedorn 
flts (solid lines), and Compound Hadron Model fits (dashed 
lines), plotted as functions of mass. 

However, the slopes in Fig. 1 are different for mesons 
and baryons, which means different Hagedorn tempera- 
tures for mesons and baryons for masses in the region 
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up to l.SGeV. This is the key observation of this paper. 
Applying the form (|^), with mg = 500MeV, we find the 
best fit for the mesons with T^^s — 197MeV, and for the 
baryons with Tbar = 141MeV. These temperatures differ 
considerably from each other, refiecting the fact that the 
baryon spectrum grows much more rapidly. 

A question arises if it were possible that meson states 
are missed in experiments much more frequently than 
baryons, such that the effect of Fig. 1 is spurious. The 
numbers show this is highly unlikely. In order to make 
the meson line parallel to the baryon line we would 
have to aggregate as much as 500 additional states up 
to TO = l.SGeV, more than the present number of 407 
states at this point. Thus, we reject the scenario of miss- 
ing meson states. We also remark that discoveries of 
some extra states can only increase the slopes in Fig. 
1, thus lowering the presently-fitted values of Tmes and 
Tbar- Another question is whether we have reached the 
asymptotics with m ^ l.SGeV. This question can only be 
addressed within a sufficiently accurate theory or model. 
It should be stressed that the function multiplying the ex- 
ponent has significance for the fitted values of Tmcs and 
Tbar- For instance, with toq = IGeV we obtain 228MeV 
and 152MeV for mesons and baryons, respectively. Using 
different values of toq for mesons and baryons may bring 
the two temperatures closer to each other. Anyway, in 
the range of to up to ^ l.SGeV we effectively have two 
distinct Hagedorn temperatures for mesons and baryons, 
as seen in Fig. 1. 
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FIG. 2. Cumulants of meson (a) and baryon (b) spectra 
with selected strangeness, plotted as functions of mass. 



In Fig. 2 we compare the cumulants of spectra for 
mesons and baryons of definite strangeness. We can see 
the universality of slopes for both mesons and baryons. 

The theoretical challenge is to explain the different 
behavior of mesonic and baryonic spectra. Historically, 
the Dual String Model [Q produced the exponentially- 
growing mass spectra. In the framework of the Veneziano 
model 1^ Huan and Weinberg |^ derived the following 
asymptotic formula: 
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where a' is the slope of the Regge trajectories, ao is the 
intercept (averaged over various trajectories), and D is 
the number of dimensions, usually assumed to be 4 
although higher values have also been attempted ||l^ . In 
the strict sense, the dual string models are constructed 
for mesons, and their application to baryons has a more 
phenomenological character. Since a' « IGeV^^ is uni- 
versal for all hadrons, the difference between mesons and 
baryons could only come from different values of ao, and 
possibly from D. We have tried all reasonable values 
for these parameters, and concluded that they result in 
very similar curves. The function multiplying the expo- 
nent in Eq. (||) has significance here. For larger D it 
decreases faster with to, which effectively compensates 
for the growth from the exponent in the interesting re- 
gion of 771. We were able to fit the meson line in Fig. 1 
with formula (jj), but it is impossible to fit the baryons. 
Moreover, the experimentally accessible range of to is 
not asymptotic enough to justify Eq. (^. Indeed, Eq. 
(^) is a direct consequence of a combinatorial problem 
known as partitio numerorum. For D = 1 Eq. is the 
asymptotic form for the number of different partitions of 
a positive integer 77. into non-negative integer components 
pT[ . For Z? > 1 it is a generalization to the case where 
the components are of D different types From the 
Regge formula (^) we can see that for to in the range 
1 — 2GeV the values of 77 lie between 1 and 4, hence 77 
is not large enough to justify the asymptotic form (jj). 
In this paper we do not pursue the dual-model track any 
further, and leave it with the conclusion that it works for 
the meson spectrum, and fails, in the present form, for 
baryons. 

Is there a simple explanation of the behavior of Fig. 
1? What we have to search for is a mechanism of the 
observed rapid increase of degrees of freedom with to. 
Inspired by the combinatorics of partitio numerorum we 
propose a candidate mechanism capable of explaining the 
result of Fig. 1 not worse than the Hagedorn model. We 
call it the Compound Hadron Model, in analogy to com- 
pound nuclei |12]. In the statistical model of nuclear 
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reactions a compound nucleus is an object whose density 
of states is enormous, and is described approximately by 
the formula p{E*) ~ f{E*) exp {h^E*^ , where E* is the 

(large) excitation energy of the nucleus, 6 is a nucleus- 
dependent constant, and f{E*) is a slowly- varying func- 
tion. Let us illustrate how this formula can arise in a 
toy example. Suppose for simplicity that we can ne- 
glect the exclusion principle, and the system has equally- 
spaced single-particle excitation levels, Sk = kAE. In the 
ground state all particles sit in the lowest level, eo = 0. 
The excitation energy can be composed in very many 
different ways. Let us express E* in units of AE, i.e. 
E* = n AE. We can make E* by exciting n times 
the lowest (first) level, E* = nei, once the second level 
and n — 2 times the first level, E* = 62 + {n — 2)ei, 
and so on. In general, E* — eiki + £2^2 + ... -I- £nfcn, 
which means n — ki + 2k2 -I- ... -I- nfc„, with integer 
ki > 0. Of course, this is nothing else but the parti- 
tion problem, and the number of distinct ways this may 
be achieved is, for a large n, given by Eq. (|) with 
D = 1. The difference is, however, that now n is pro- 
portional to the excitation energy, and not the square of 
mass, as in Eq. (^). One can extend this combinatorial 
problem to account for D different types of excitations 



/c|"^ > 0, which gives the number of partitions, Poin), of 
Eq. (|^). In addition, we may incorporate the exclusion 
principle by disallowing the repetition of the same level, 
i.e. requesting that none of the numbers /c^^"'' is repeated. 
This leads to the asymptotic formula [fsl 
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where "NR" indicates "no repetitions" . Note the factor 
of under the square root compared to ^ in Eq. (|^) 
(there are less partitions without repetitions than with 
repetitions), and the absence of D in the power n multi- 
plying the exponent. 

Inspired by the compound nucleus model and by the 
expression for P^^(n), we propose the following guess 
formula for the hadronic mass spectra, which is a direct 
consequence of Eq. (||): 
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where c is a constant, ttiq is the ground-state mass, D is 
the number of types of excitations, and AE' is the aver- 
age level spacing. The asymptotic power of the function 
multiplying the exponent is —5/4, whereas the power 
—2 is used in an analogous formula in compound nu- 
clei [12[. The underlying physical picture is following: 



hadrons are bound objects of constituents (quarks, glu- 
ons, pions). The Fock space contains a ground state, and 
excitations on top of it. In the case of the compound nu- 
cleus these elementary excitations are Ipl/i, 2p2h, 3p3/i, 
etc. states. In the case of hadrons they are formed of 
qq and gluon excitations, e.g. for mesons we have qq, 
qqg, qqqq, qqgg, etc. Now, we form the excitation energy 
(hadron mass) by differently composing elementary exci- 
tations, which bring us to the above-described partition 
problem. It seems reasonable to take zero ground-state 
energy for mesons, m^'^^ = 0, which are excitations on 
top of the vacuum. For baryons we take mg'^'' = 900MeV, 
which is basically the mass of the nucleon. The quantity 
AE/D is treated as a model parameter and is fitted to 
data. 

The results of the compound- hadron- model fit, Eq. 
(^), are shown with the dashed line in Fig. I. The curves 
are slightly bent down, compared to the Hagedorn-like 
fits (solid lines), which is caused by the square root in 
the exponent of Eq. (^). But the fits are no worse, or 
even better when the fit region is extended to m = 2GeV. 
Numerically, the least-square fit for m up to l.SGeV gives 
AE'^'^yO = 50MeV for mesons, and AE^'^'/D = 53MeV 
for baryons. The proximity of these numbers shows that 
the scales for mesons and baryons are similar. The num- 
ber of types of the excitations, D, is taken to be I. This 
means that we treat all levels, carrying spin, orbital, ra- 
dial, or isospin quantum numbers, as separate. If we 
treat them as D-times degenerate, then the average spac- 
ing AE increases D times, such that the dependence 
on D cancels. Note that in case of gluonic or qq exci- 
tations color does not bring extra degeneracy, since all 
states have to be color-neutral and there is only one way 
in which the elementary excitations can be coupled to 
a color singlet. The obtained values for Ai?™°*' (with 
D = 1) mean that the corresponding n aXm — I.SGeV is 
around 35 for mesons and 17 for baryons. Such n are suf- 
ficiently large to justify the use of the asymptotic formula 

There is an interesting effect we wish to point out. It 
is natural to expect that a bound hadronic system has an 
upper limit for the excitation energy. It is helpful to think 
here of bags with a finite depth, or of breaking-up flux- 
tubes. Thus, in constructing the Fock space for bound 
objects we should have a limited number of quanta to our 
disposal. Translating this idea into our toy combinatorial 
problem, we now have to ask about the number of ways 
of partitioning n (without repetitions) into integer com- 
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now with /cmax > fci"'' > 0, where /cmax is the maximum 
available number. We denote this number of partitions 
with limited k'f^ as P^^(n, /cmax)- Figure 3 shows the 
decimal log of cumulants of this quantity for D — \, 
^"^j^ Pj^^(i, fcmax), plotted as a function of n for vari- 
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partitions with unlimited kl"'. Each curve becomes flat 
at n = A;max(fcmax + l)/2, which is the point where all 
available levels have been filled. Beyond that point we 
cannot form bound states any more. We term this effect 
combinatorial saturation. It means that in a compound- 
hadron model beyond a certain mass m no more light- 
quark resonances can be formed. For instance, our toy 
model with D = 1 and fcmax = 12 predicts that this hap- 
pens at n = 78. With AE obtained in the fits, this leads, 
according to the formula m = mo + nAE, to a maximum 
mass of ~ 4GeV for light-quark mesons and ~ 5GeV for 
light-quark baryons. We can see in Fig. 3 that before 
reaching the combinatorial saturation, the curves with 
a given /cmax depart from the unlimited-partition curve, 
^max = cxD. It turns out that the fcmax — oo curve is very 
close to the asymptotic (large n) curve, described by for- 
mula (^, even for low values of n. This is why Eq. (||) 
can be used even at low values of m. We may thus spec- 
ulate that the departure of data from the dashed lines in 
Fig. 1 around m = 2GeV, similar to the behavior of Fig. 
3, is a signal of the combinatorial saturation. More com- 
plete data in that region and above would definitely help 
to check if this is really the case. We note that Bisudova, 
Burakovsky and Goldman fl^ obtain similar limits for 
mass of light hadrons from non-linear Regge trajectories. 
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FIG. 3. Cumulants of Pf^(i,fcn,ax) for various values of 
fcmax- The integers at the right ends of the curves are the 
saturation values of the cumulants. 

Let us summarize the basics of combinatorics behind 
the formation of the hadronic spectra. As we increase the 
mass, more elementary excitations may participate in the 
resonance-formation process, in other words more degrees 
of freedom open up. This, via partitio numerorum, leads 
to an exponential growth of the number possibilities of 
forming excited states. Beyond a certain mass there are 
no more elementary excitations available, and we reach 
the combinatorial saturation in the number of hadronic 
states. Clearly, the real life of hadronic physics is more 
complicated than the toy model discussed above. How- 
ever, the result may be robust and not care about the de- 
tails, for example whether the levels are equally spaced. 



as assumed here, or not. The success of the statisti- 
cal model in nuclear physics shows that details are not 
needed to understand gross features system with many 
degrees of freedom. 

Finally, let us comment on the different power of m in 
the exponent of Eq. (||) and the Hagedorn hypothesis, 
Eq. (^. It is not a paradox, since our simple assumptions 
in the language of quark and gluon (or pion) excitations, 
are completely different than the assumptions of the Sta- 
tistical Bootstrap Model, where self-similarity of fireballs 
is the key ingredient. The thermodynamical implications 
of the two models are different: whereas in the Hagedorn 
model the temperature T of a hadronic system cannot 
be larger than the Hagedorn temperature, in the Com- 
pound Hadron Model there is no formal limit on T. We 
stress that this is of formal significance only, since de- 
confinement is expected to occur at temperatures of the 
order of ^ 150MeV, dissolving hadron resonances into 
the quark-gluon plasma. 

One of us (WB) is grateful to A. Bialas, A. Horzela, J. 
Kwieciiiski, and P. Zenczykowski for useful discussions. 
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